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Abstract 

We  consider  a  space-time  variational  formulation  for  linear  parabolic  partial  differential  equations.  We  introduce 
an  associated  Petrov-Galerkin  truth  finite  element  discretization  with  favorable  discrete  inf-sup  constant  j3g\  /3s  is 
unity  for  the  heat  equation;  (3g  grows  only  linearly  in  time  for  non-coercive  (but  asymptotically  stable)  convection 
operators.  The  latter  in  turn  permits  effective  long-time  a  posteriori  error  bounds  for  reduced  basis  approximations, 
in  sharp  contrast  to  classical  (pessimistic)  exponentially  growing  energy  estimates. 

Resume 

Nous  considerons  une  formulation  variationnelle  espace-temps  pour  les  equations  differenticlles  paraboliques 
lineaires.  Nous  y  associons  une  discretisation  par  elements  finis  de  Petrov-Galerkin  pour  laquelle  la  constante 
de  stability  inf-sup  (3g  possede  des  proprietes  agreables  :  (3s  est  unite  pour  l’equation  de  la  c.haleur;  f3g  a  une 
croissance  seulement  lineaire  en  temps  pour  des  operateurs  de  convection  non-coercifs  (mais  asymptotiquement 
stables) .  Dans  le  cadre  des  approximations  par  bases  reduites,  cette  derniere  propriete  permet  d’obtenir  des  bornes 
efRcaces  pour  l’erreur  a  posteriori  en  temps  long,  en  net  contraste  avec  les  estimateurs  d’erreur  en  energie  classiques 
(pessimistes)  qui  presentent  une  croissance  exponentielle. 
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Version  franchise  abregee 

Soit  l’equation  aux  derivees  partielles  parabolique  (1)  et  sa  formulation  variationnelle  espace-temps  : 
u  G  X  :=  L2{I  =  (0,T];U)  D  H^I-V')  :=  {v  G  :  v(0)  =  0},  verifie  (2)  pour 

l’espace  des  fonctions  test  y  :=  L2(I;V);  ici  b( w,v)  :=  JjKwit),  v(t))v'xv  +  a(w(t),v(t))]dt  et  f(v)  := 
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fj(g(t),v(t))v'xvdt.  Les  normes  sur  X  et  y  sont  definies  par  \\w\\2x  :=  ||u>  IIL(W+II-IIL(^)+II-(t)IIh 
et  |M|y  :=  ||^||l2(/;V)-  On  considere  V  =  Hq(SI)  (operateurs  spatiaux  de  2e  ordre),  ou  f 2  est  le  domaine 
spatial.  On  montre  (Proposition  1  basee  sur  [6])  que  pour  des  problemes  coercifs,  la  constante  inf-sup  0 
de  b  est  positive  et  bornee  inferieurement. 

Nous  introduisons  ensuite  une  discretisation  par  elements  finis  de  Petrov-Galerkin  :  soit  us  G  Xs 
verifiant  b{us,vs)  =  f{vs),\/vs  G  y, 5;  avec  Xs  :=  SAt  <8>  Vh  et  ys  =  QAt  '®  Vh,  S  =  (A t,h),  ou  SAt,Vh 
et  Qm  sont  respectivement  des  elements  finis  lineaire  par  morceaux  en  temps  (de  pas  At)  et  en  espace 
(de  diametre  h),  et  constants  en  temps  par  morceaux.  Cette  methode  coincide  avec  une  discretisation 
de  Crank-Nicolson  en  temps,  ainsi  nos  resultats  s’appliquent  directement  a  ce  schema  standard.  Nous 
definissons  une  norme  modifiee  sur  A’,  |||  •  |||ar,,5,  ou  la  partie  fx  ||in||ydt  de  |  •Ik  est  remplacee  par 
une  somme  des  moyennes  de  w  sur  chaque  element  temporel.  Nous  montrons  (Proposition  3)  que  pour 
l’equation  de  la  chaleur,  les  constantes  inf-sup  discrete  0s  et  continuite  discrete  7 5  sont  unites. 

Nous  considerons  maintenant  une  approximation  par  bases  reduites  un,  comine  dans  [2].  Nous  pouvons 
construire  a  la  Proposition  5  des  bornes  sur  l’erreur  a  posteriori  pour  le  champ  solution  un  et  pour  la 
sortie  scalaire  Sjy.  L’utilite  de  ces  bornes  est  tres  liee  a  la  dependance  de  0s  aux  parametres  du  probleme 
et  au  temps  final  T.  La  procedure  de  calcul  de  ces  bornes  est  similaire  a  la  decomposition  hors  ligne-en 
ligne  des  bases  reduites  standard. 

Enfin,  nous  presentons  des  resultats  numeriques  pour  la  constante  inf-sup  discrete  0s-  Pour  des  opera¬ 
teurs  non- coercifs  rnais  asymptotiquement  stables  (e.g.  convection),  on  observe  que  le  parametre  inf-sup 
decroit  seulement  en  (/iqT)-1,  ou  /j-|  est  la  vitesse  de  deplacement  et  T  le  temps  final.  L’exemple  est  un 
domaine  spatial  unidimensionnel  (0,1)  avec  l’operateur  p\(x  —  \/2)ux  —  uxx  pour  pi  >  2ir2  (/xi  <  27r2 
entraine  un  probleme  coercif).  Les  estimateurs  associes  sont  beaucoup  plus  precis  que  ceux  classiques  en 
energie  qui  predisent  une  croissance  exponentielle  en  e^lT  et  qui  sont  done  inutilisables  en  pratique. 


1.  Space-time  formulation 

We  first  formulate  a  general  linear  parabolic  equation.  Consider  Hilbert  spaces  V  ^  H  ^  V'  and 
an  operator  A  G  jC(V,V'),  (A<j>,ip)v>xv  =  f°r  W  G  L  Setting  /  :=  (0,T],  T  >  0  and  given 

g  G  L2(/;  W),  we  look  for  u  such  that 

u(t)  +  Au(t)  =  g(t )  mV',  t  G  I;  u( 0)  =  0,  (1) 

and  an  associated  output  of  the  form  s  :=  fj£(u(t))dt  for  some  i  G  V' .  We  restrict  ourselves  to  LTI 
systems  even  though  some  of  our  results  can  be  extended  to  a  more  general  situation. 

For  the  derivation  of  the  space-time  variational  form  of  (1),  we  introduce  the  trial  space  X  := 
L2(I;V)r\H{0)(I-,V'),  where  H{0){I;V')  :=  {u  G  V7)  :  v(0)  =  0}  with  norm  \\w\\%  :=  \\w\\l2(I.v)  + 
\\w\\l2{i.v>)  +  llu;(^1)llff  and  the  test  space  y  :=  L2{I\V)  with  norm  Ill’ll^  :=  |M|z,2(/;iq.  Note  that 
X  =  {L2{I)  ®b)n  ®  V ')  and  y  =  L2(I)  ®V  which  will  allow  for  a  tensor  product  discretization. 

Then  definitions  b(w,v)  :=  Jj[(w(t),  v(t))vrxv  +  a(w(t),v(t))]dt  and  f(v)  :=  fj(g(t),v(t))v'xvdt  yield 
the  space-time  variational  formulation 

u  G  X  :  b(u,v)  =  f(v),  \/v  G  y.  (2) 

The  well-posedness  of  (2)  has  been  shown  (under  suitable  assumptions)  in  [6,  Theorem  5.1]. 

The  approach  of  [6]  can  also  yield  an  estimate  for  the  inf-sup  constant  0  :=  inf^g^  sup„g-y  y  ■ 

We  define  g  :=  sup0^ey  |||^  and  0*  :=  inf^gy  sup^eV  w^HWv  ;  we  then  claim 
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Proposition  1  Assume  that  there  exist  Ma  <  oo,  a  >  0,  and  A  >  0  such  that  \a(tp,  4>)\  <  Ma  ||4||v  W^Wv 

(continuity)  and  +  A||^||^-  >  a||^||y  (Carding)  for  v,w  £  V.  Then,  we  obtain  the  inf-sup  lower 

hound  ft  >  flLB  —  min{l,(g-Ag-)min{l,A/-2}} 
uounap^p  —  max{i,(/3*)— 1}v/2 


Proof.  (Sketch)  Let  0  ^  w  £  X.  Set  zw  :=  (A*)  1w  (where  A*:  V  — >  V'  denotes  the  adjoint  of  V)  and 
set  vw  :=  zw  +  w  £  y.  Then,  <  2 max{l,  (/3*)-2}||u>|||,.  Using  the  estimates  ||u>(f)||y/  < 

Ma\\zw(t)\\v,  (w(t),zw(t))vxv  >  (a  —  Ap2)M“2||i(;(<)||y,  and  a(w(t),  zw(t))  =  \ft\\w(t) \\2H  we  arrive  at 
b(w,vw)  >  (a-A02)(M-2||w||22(/;y,)  +  ||w||22(/;V))  +  ||«;(T)||2f  >  mm{(a- Xg2)mm{l,  M~2},l}\\w\\2x  > 

/3LB||u>||,r  |K||y,  which  proves  the  claim.  □ 

Remark  2  Note  that  (3LB  does  not  depend  on  the  final  time.  However,  the  estimate  is  only  meaningfid 
if  a  >  A g2 ,  i.e.,  if  the  system  is  coercive.  In  the  non-coercive  case,  (1)  is  often  transformed  via  u(t)  := 
e~xtu(f)  to  obtain  a  coercive  problem;  however,  this  yields  an  inf-sup  bound  that  behaves  as  e~XT  —  often 
extremely  pessimistic  and  clearly  unsuitable  for  error  estimation  in  long-time  integration. 


2.  Petrov- Galerkin  Truth  Approximation 

Let  Xg  C  X,  y$  C  y  be  finite  dimensional  subspaces  and  u$  £  Xg  the  discrete  approximation  of  (2), 
i.e.,  b(ug,vg)  =  f(vg),  Vv$  £  y$,  sg  =  /QT  £(ug(t))  dt.  Henceforth,  we  concentrate  on  the  case  H  =  L2(U), 
V  =  Hq  (f2).  Let  Xs  =  5a*  0 14,  34  =  Qa*  0  14,  S  =  (At,  h),  where  5a*,  14  are  piecewise  linear  and  Qa* 
piecewise  constant  finite  elements  with  respect  to  triangulations  7)fpace  in  space  and  Ta™6  —  {U-1  = 
(i  —  1)  At  <  t  <  i  At  =  tl,  1  <  i  <  r}  in  time  for  At  :=T/r.  This  coincides  with  the  Crank-Nicolson  (CN) 
scheme  if  a  trapezoidal  approximation  of  the  right-hand  side  temporal  integration  is  used;  hence,  we  can 
derive  error  bounds  for  the  CN  scheme  via  our  space-time  formulation. 

We  introduce  a  different  norm  on  X:  For  w  £  X  and  1 1  :=  (tl~l,tl),  set  wl  :=  (At)-1  ffi  w(t)  dt  £  V 
and  w  :=  XP  0  do1  £  L2(/;  U),  where  \P  is  the  characteristic  function  on  P;  then,  set  |||«;|||^.  s  := 

IHIi2(/;V')  +  IMI|2(i;v)  +  \\w(t)\\h  and  the  inf-sup  parameter  /3g  ■=  '^w5exs  sup,,^^  111^1^1^10  as 
well  as  the  stability  parameter  :=  sup.^^  supt,je:yj  iii^qj^lliiy  • 

Proposition  3  Let  a(-,-)  be  symmetric,  bounded  and  coercive  and  set  \\4>\(y  :=  a(4>,(j>),  f>  £  V;  then  we 
have  (3g  =  7<5  =  1  • 

Proof.  (Sketch)  Since  vs  £  ys  is  piecewise  constant  in  time,  we  have  Jj  a(ws,vs)  dt  =  a(uig,vg)  dt  for 
all  ws  £  Xs-  Hence,  b(wg,vs)  =  jIa(A~h1wg  +  wg,vg)dt,  where  zg  :=  Af^wg  is  defined  by  a(zg,<f>h)  — 
( wg,(f>h}v'xv >  £  Vh-  Note  that  =  || || v'' -  We  may  prove  that  for  any  vs  £  ys  there 

exists  a  unique  zs  £  Xg  such  that  f/a(Af1zg  +  zg,qg)dt  =  fz  a(vg,  qs)  dt  for  all  qs  £  34-  Note  that 
vg  :=  Af^zg  +  zg  £  yg  for  zg  £  Xg  and  we  obtain  b(wg,  vg)  =  fx  a(Af1wg  +  wg,  Aflzs  +  zg)  dt  so  that 
sup„je  irar  =  suPvseys(Jia(Ah1™s+™fi’Ah1zs+zs)dt)/\\A((1zg+zg\\y  =  WA^wg+uisWy.  The  fact 
that  ||  Ah  =  ||  Ah  l2(iv)  ’d,s)v,xv  dt.  =  ||  wg  ||  £2(/.yq  "U 1 1  ^5 1 1  z>2  (_r- v^)  "t" 

INCOI&HIKII&.  s  proves  the  assertion.  □ 

Remark  4  Proposition  3  also  shows  the  well-posedness  of  the  discrete  problem  with  continuity  and  inf- 
sup  constant  being  unity. 
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3.  The  Reduced  Basis  Method  (RBM) 


Now,  let  p  =  {pi,p2)  G  V  :=  K2  be  a  parameter  vector  and  A  =  A(p)  :=  —Au+p\  (3(x)-\/u+p2  u,  i.e., 
a  diffusion-convection-reaction  operator  with  convection  field  (3.  We  then  introduce  a  standard  Reduced 
Basis  (RB)  approximation  [1,4,5]  for  the  Crank-Nicolson  interpretation  of  our  discrete  problem,  un{p)  G 
%At,N  =  SAt  ®  Vat;  here  Vn  C  Vi,  is  an  RB  space  provided  for  example  by  the  POD-Greedy  procedure  of 
[2].  The  RB  output  is  then  given  by  sjv  =  fT  l(v,N(t))dt(=  Jj£(uN(t))dt).  (It  is  possible,  alternatively,  to 
consider  a  space-time  RB  approximation  as  well  [7].)  It  is  then  simple  [5]  to  demonstrate 
Proposition  5  The  RB  error  satisfies  |||w,5  —  w/v|||;t,5  <  llrw \\y' / fiY B>  where  rjv(w)  :=  f(v)  —  b(uN,v), 
Vu  G  y,  and  figB  is  a  lower  bound  for  the  fig  defined  in  Proposition  3.  Furthermore,  |s<s  —  sjv|  < 

Mv'VT\\rN\\y'/rtB. 

The  utility  of  these  a  posteriori  error  bounds  is  critically  dependent  on  the  dependence  of  fig  as  a  function 
of  the  parameter  p  and  final  time  T,  fig(p;T). 

We  briefly  comment  on  the  computational  implications  of  this  space-time  error  bound.  As  always,  our 
model  problem  (as  is  the  case  here)  must  honor  the  “affine-in- functions-of- parameter”  condition  [5]  to 
permit  effective  offline-online  computation  of  the  reduced  basis  approximation  and  a  posteriori  error 
bounds.  The  current  formulation  introduces  the  further  complication  of  the  space-time  norms.  In  fact, 
calculation  of  the  dual  norm  ||  •  ||;y/  (through  the  corresponding  Riesz  representation  [5])  does  not  couple 
different  time  steps/elements,  and  hence  the  additional  online  difficulty  is  relatively  minor.  (The  offline 
stage  for  computation  of  the  inf-sup  lower  bound  by  the  Successive  Constraint  Method  [5,8]  does  require 
special  treatment,  in  particular  due  to  the  norm  |||  •  |||a’,(5-) 


4.  Numerical  Results 

We  report  numerical  results  for  the  Crank-Nicolson  scheme  for  various  choices  of  the  parameters  / j.  i ,  p2 
as  well  as  for  different  time  steps  At  and  uniform  mesh  sizes  h.  For  simplicity,  we  consider  the  univariate 
case  (in  space)  =  (0,1)  and  choose  (3(x)  =  x—\.  Let  us  denote  by  fig(p;T),  7 g(p;T)  the  numerical 
values  for  the  truth  inf-sup  and  continuity  constants,  respectively,  corresponding  to  parameter  p  and  final 
time  T.  We  start  by  confirming  Proposition  3.  Thus,  we  choose  pi  =  p2  =  0;  for  several  values  of  T,  h, 
and  At  we  invariantly  obtain  1.000  for  both  fig(p;T)  and  7 g(p-,T),  as  must  be  the  case. 

Next,  we  investigate  the  case  of  convection,  p2  =  0,  in  which  case  a  is  coercive  only  for  p\  <  2tt2. 
We  obtain  fis({ 50,0);  1.0)  =  0.050,  ft  ((50, 0);  2.0)  =  0.027  (75  =  2.60)  and  ft((100, 0);  1.0)  =  0.019, 
ft((100, 0);  2.0)  =  0.010  (7 g  =  5.6).  These  numbers  are  relatively  invariant  for  sufficiently  small  h  and  At. 
We  observe  numerically  an  overall  behavior  of  fig((pi,  0);  T)  ~  (piT)-1  and  75 ( (/^i ,  0);  T)  ~  pi  (the  latter 
is  readily  proven,  but  not  the  former).  Note  T  =  0(1)  is  effectively  a  “long  time”  in  convective  units, 
1/pi-  We  emphasize  that  although  the  problem  is  non-coercive,  the  problem  is  asymptotically  stable 
in  the  sense  that  all  eigenvalues  er  of  —a(i[>,<l>)  =  er(^>,  <t>)v'x.v  lie  in  the  left-hand  plane;  this  stability 
is  reflected  in  the  inf-sup  behavior.  In  contrast,  a  standard  energy  approach  [3]  gives  effective  inf-sup 
constants  on  the  order  of  e-MlT  (here  about  10-8).  Hence,  the  traditional  method  fails  to  provide  useful 
results,  whereas  our  new  approach,  which  reflects  the  true  time-coupled  properties  of  the  system,  yields 
relatively  sharp  error  bounds. 

Finally,  we  consider  the  case  p\  =  0  which  gives  rise  to  an  asymptotically  unstable  (and  non-coercive) 
system  for  p2  <  — 7r2.  This  means  that  any  error  estimate  must  grow  exponentially  with  the  final  time 
T.  We  observe  this  for  our  estimator  as  well,  e.g.  ft((0,  —20);  1.0)  =  3.93  •  10-5  (order  of  e^2T). 
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